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الباب الاول 


مقدمة في التحلیل العددي 
1 -1 استهلال introduction‏ 


التسمية الدقيقة لهذا ااتخصص هو التحلیل العددی( numerical‏ 
(analyses‏ بستخدم التحلیل العددي في حل المعادلات الرياضية التي يصعب حلها 
أو يستلزم وقت طویلا في الحل. فإذا كانت لدینا معادلة رياضية من الدرجة 
الرابعة أو الخامسة فان من الممکن حلها بالطرق التقليدية. و حتی هذه الطرق 
ناخد وقتاء اما 13 كانت المعاددة من الدرجة السادسة و آکثر یکون من الصعب 
Laila‏ التعامل مع المعادلة تقلیدیاء و مع ظهور الحاسب الآلي ( الکومبیوتر) 
اتضحت آهمیته البالغة في حل هكذا معادلات. و ذلك لتوفیره الوقت و الجهد. و 
بالاخص في المعادلات التي تحتاج إلى تکرار كبير من Jal‏ الوصول إلى النتيجة 
أو الحل, 


عندما یحصل فني أو مهندس على وظيفة معینة Mia‏ في الصيانة أو 
الإنتاج أو الجودة أو تدبير al gall‏ الأولية» يتوقع أنه سيفتح جهاز الحاسوب لیجد 
iLi‏ أحد برامج SAP‏ أو Baan‏ أو si Oracle si PeopleSoft‏ 
JDEdwards‏ و کل ما عليه هو الضغط على بضعة آزرار هنا و هناك لیعرف 
dia‏ المخزون و الطلبیات و موعد تسلیمها و نتانج اختبارات الجودة و كمية 
المواد الأولية التي يجب شراو‌ها و عدد ال lots‏ في خطوط الانتاج و أسماء 
العمال الذین عملوا في تجمیع منتج معین و عدد ساعات غیاب عامل ما و 
الاقتراحات التي قدمها العمال و عدد توقفات АЙ‏ ما ...لکن هذا غير موجود على 
آرض الواقع. واقعنا العربي... 


و حتى الشركات الغنية كشركات البترول و الصلب. التي استثمرت في هذه 
الحلول و ركبتهاء لا تستفيد منها كاملا و أغلب الموظفين يعتبرونها عبئا إضافيا 
على واجباتهم اليومية و قد يستعملونها فقط عند قدوم زوار خارجيين. أما في 
الشركات المتوسطة و الصغيرة. فمن الصعب على الفني أو المهندس إقناع 
المدراء بالفائدة الملموسة لمثل هذه الأدوات خصوصا أن ثمن تملكها السنوي 
Total Cost of Ownership‏ مرتفع مقارنة بميزانية الشركة السنوية لذلك 


فالحل الوحيد الذي یجده المهندس لتدبیر آموره اليومية هو ابتکار حلول برمجية 


و لن يتأتى له ذلك الا بمعرفة مسبق بطرق التحلیل العددي و استعمال 
الحاسوب في حل معادلات رياضية. و هذا يدرس في الجامعات و مراکز التکوین 
з)‏ منتدانا | لحبيب) و ليس في مكان العمل. OY‏ وقت العمل ضيق و محسوب.كما 
آنه لا مجال , لاستعمال برامج dive Дл‏ تحت طائلة تعریض الشركة لغرامات مالية 
في حال قيام مايكروسوفت و أخواتها بالتحقق من برمجيات الشركة fil)‏ مكالمة 
هاتفية من عامل لم يستفد من العلاوة السنوية ). 


ور Eb eus ee‏ و وی میور 
تستعملها à‏ لمنزل لا تس (xn‏ با كلها في حاسوب العمل دون 














في النھایف آقول اني لا آقصد أن یکون المهندس مبرمجا محترفا زيادة 
على التخصص الذي درسه. لکن وضع برنامج صغیر ل نبع مؤشر ما او ورقة 
AI‏ مب نت الرياضية أسرع و آنجع من تعلم برامج من العیار الکبی 











في حقيقة الأمر یوجد الکثیر من البرامج یمکن استخدامها في هذا المجال 

مثل برنامج MATHCAD , MATHLAB‏ و حتى برنامج EXCEL‏ 
NUM pue‏ 097 يەكن أجراءه في حل الكثير من سای 3 ч‏ 
لفات ا في التعامل مع المع ادلات الرياضية اس Lu га PA‏ 
Lal FORTRAN‏ تمناز به هذه ARTI‏ من استقرار عال و دقة في المعادلات و 
تحدید دقیق لنوع المتغیرات مع سعة كبيرة في نوع المتغیرات» من الممکن آنها 















ليست على Аз уз‏ كبيرة من الانتشار و الاستخدام + غير أن المتخصصین في مجال 
الت ل العددي و الریاضیات بصفة عامة یعترفون لها بالفضل نستخدم ЛУ!‏ 





الرابع منهاو هي 
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و هو إصدار يعمل تحت بيئة الوبندوز WINDOWS)‏ ) و لست في معرض 
الدعاية لهذه اللغة و لا أحب المقارنة بين اللغات المختلفة غير أن 


FORTRAN LANGUAGE‏ هي لغة المهندسين بكل كفاءة. 


пани сна ia ika | ہو بعض‎ beum غير أن هذا‎ 











قبل الت ET‏ المعدلات الرياضية 
التواب و المتغيرات. 


أساسيات لغة فورترا FORTRANO‏ 








CHARACTERS الرموز‎ 1 


- 





- 


NUMBERS CHARACTERS الأرقام‎ 1 
9 87 6 5 4 3 21 0 و هي‎ 
ALPHABATIC CHARACTERS الحروف‎ 1.2 











21 آنواع البیانات في 44 FORTRAN‏ 


لغة الفورتران FORTRAN‏ باحدی الأسالیب الاتية : 








INTEGAR صحيحة‎ 















طرق حل العادلات ذات المجهول الوا 





JATION IN 








l-2‏ طريقة المقاطع 
BISECTION METHOD‏ 
لنفرض أن Gal‏ المعالة الرياضية التالية : 
0 - 4 + 2۸2 


مع ملاحظة أن الرمز(۸) يعني أس أي مربع (X ) АА‏ و 2۸45 يعني الأس 


سيكون الحل المثالي لهذه المعادلة هو كالأتي: 
X^2-4‏ 
Хан‏ 
X-2‏ 
за sX=-2‏ حل دقیق تماما بمجرد إتباع الخطوات السابقة, نحصل على الحل. 
و إذا كانت لدينا المعادلة الرياضية التالية : 
Х^3+4Х^2-10=0‏ 


فحلها يرهق قليلا و لا يعطي قیم دقيقة » و لهذا نلجأ إلى طريقة المقاطع 
BISECTION METHOD‏ التي تتلخص Led‏ يأتي : 


لنفرض أن لدینا F(X) Alla‏ معرفة و مستمرة خلال الفترة } A, В‏ { مع 
الدالتین F(B) , F(A)‏ فهنالك قيمة ХО‏ تنتمي إلى الفترة 


F(X) التي تکون حلا للمعادلة‎ 5: F(X0)=0.0 Walls {A,B} 
تنتمي إلى نفس الفترة فهذا يعني ان الدالتین‎ А В بما أن القيمتين‎ 


F(B) , F(A)‏ تنتمي إلى نفس النطاق و يكون الرسم البياني للدالة کالاتي: 








FIONT SOLUTION 


(شکل 2 -1) 
و حل المعادلة یکون هو النقطة PIONT SOLUTION) ХО‏ ) و ہما أن 
آساس الحل تقريبي فتکون الطريقة المثالية لخطوات dall‏ هي اخذ مجموعة من 
النقاط التقريبية المقربة للقيمة الصفرية كالقيمة 0.00001 
و اختبارها فهي لا تساوي صفر لکنها دقيقة جدا و مقربة جدا من الصفر حيث 
يصعب الحصول على الحل الصفري في التحليل العددي مباشرة و هذه القيمة 
سنفرضها ليحملها المتغير P‏ و يقترب كل مرة بقيمة دقيقة من الحل الصفري و 
تعتمد طريقة المقاطع (BISECTION METHOD)‏ على نفس الفكرة في 
التعامل مع المسائل الرياضية و خوارزمية البرنامج هي كالآتي : 
START‏ 
RAED(A,B)‏ 
DO P=(A +B )/2‏ 
Y1= A^3 + 4*А^2 - 0‏ 
Y2- ВАЗ + 4*B^2 10‏ 


У0=Р^3 + 4*P^2 – 10 


IF (Y0*Y1> 0.0) GOTO STEP 10 





DO A=-P 
GOTO STEP 20 








SOLUTION C THIS PROGRAM TO CALCOLATE 
METHOD OF EQUATION BY USED BISECTION 














RAED (*,*) A,B 
P=(A +В 2 

У1=А^3 + 4*А^2 – 10 

У0=Р^3 + 4*P^2 – 10 





А=А 
В=В 
А=-Р 





ЕТТЕ(*,*) У0,Р 


STOP 











1 بها القيمة‎ = en هر لد لك‎ à: 





Nuten method طريقة نيوتن‎ 2-2 


في هذا الدرس نستعرض سويا احد الطرق المتبعة في التحلیل العددي بعد 
أن تکلمنا عن مفهوم التحلیل العددي و آهمیته و الطريقة الثانية المتبعة في 
التحليلي العددي هي طريقة نیوتن newten method)‏ ( و لا آود أن أتطرق 
سی الاستتتاج الرواطسي لهذه الطريقة بالق الذي ارغب فيه في الترکیز على 
البرامج الرياضية لها, تعتبر هذه الطريقة من الطرق السهلة في إيجاد القيم 
التقريبية للمعدلات الرياضية و الاختلاف الأساسي بيناه و بين الطريقة السابقة ( 


طريقة المقاطع) (bisection method)‏ هو سهولة استخدام طريقة نيوتن و 
سرعتها في إيجاد الحل التقريبي للمعادلة. 


المفهوم الرياضي لنيوتن 


اذا كانت لدينا Alla‏ حقيقية f(x)‏ و كان 0× هو الجدر للمعادلة المطلوب الحصول 
عليه ( الحل التقريبي للمعادلة) فيمكن إيجاده عن طريق الآتي: 

X2=x1*f(x)/f '(x) 

حيث : 

2 هي القيمة التي نبحث عنها أي نفس قيمة الجدر ( قيمة (х0‏ 

1 هي القيمة المدخلة عند القراء 

F(x)‏ هي قيمة المعادلة بعد تعويض بي قيمة x]‏ في المعادلة 

F ' (x)‏ هي قيمة المشتقة الأولى للمعادلة بعد. 

و تكتب خوارزمية البرنامج كالأتي: 


Read хі 1 
#(х1)=х1^3+4*х1^2-10 2 


f ' (x1) = 3*x1^2 +8*х1 3 

if f(x1) > 0.00001 write x2,f(x1) 4 

x2=x1-( f(x1)/ f ' (x1) ,xl=x2 goto .5 
step (2) 

write х2 , .6 

stop .7 


شرح الخوارزمية 

يبدأ البرنامج بقراءة قيمة الدالة عند النقطة x1)‏ ) ثم تتم عملية تعويض من dal‏ 
إيجاد القيمة الفعلية للدالة عند نفس النقطةء تعوض قيمة х1)‏ ) في مشتقة الدالة و 
كأننا أوجدنا الميل في هذه АЛАЛ‏ تتم اختبار قيمة الدالة f (х1)‏ فاذا كانت أصغر 
من 0.00001 فتتم طباعة القيمة مباشرة و إلا ننتج قيمة جديدة هي 2× من 
طرح قيمة 1× من مقسوم الدالتين sf'(x) ,f (x)‏ هي قيمة أو نقطة تقاطع 
المماس مع محور السينات و عندما نجعل قيمة x] =x2‏ فإننا نقترب من الحل JIS)‏ 
أي من النقطة الصفرية التي تحقق الحل. 


و هذا هو البرنامج بلغة --FORTRAN‏ 


Read(*,*) ×1 

10 f(x1)=x143+4*x1%2-10 
f ۲ (х1) = 3*x1^2 +8*х1 
if f(x1) > 0.00001) goto 20 
x2-x1-(f(x)/f ۲ (x)) 
xl-x2 

goto 10 

20 Write(*,*) x2 р 
stop 
end 


و یمکن استخدام рО La‏ بدلا من جملة GOTO‏ لأداء نفس المهمة على 
النحو الآتي: 


Read(*,*) x1 e 
DO 500,I=1,N • 
#(х1)=х1^3+4%*х1^2-10 ° 


f'(xi)- 3*x1^2 +8*х1 e 

if f(x1) > 0.00001) GOTO 20 e 
x2=x1-(f(x)/f' (x)) є 

xl-x2 e 

GOTO 500 • 

500 CONTINUE • 
20 Write(*,*) x2,p 9 
stop ° 

end e 


و الرسم التالي یوضح فكرة نیوتن في حل المعدلات : 





= 
= 


و نلاحظ انه بدأ من نقطة بعيدة و في کل مرة يقترب بمقدار معین من القيمة 
المطلوبة X)‏ و كأنه يأخذ الشکل التالي: 


الهدف( لقيمة الصفریة) 


d 





(شکل 3-2( 


2 -3 طريقة نیوتن رافسن ونیوتن رافسن المعدلة 


Newton-Raphson Method 
And 
Modified Newton-Raphson Method 


طریقتا نیوتن رافسن ونیوتن رافسن المعدلة: 
هما طریقتان لحل المعادلات الانية غير الخطية , ولو أخذنا معادلتین في 
متغیرین كمثال ونرید ایجاد قيمتي х,у‏ اللتان تحققان المعادلة )1( 
F1(x,y)=0 F2(x,y)-0............ (1)‏ 

آولا: طريقة نیوتن رافسن 
من متسلسلة تایلور في сы мала‏ کمثال 

Еу) EEL yA С AGO سل‎ 

172 50.7777200 E (кыек) کرک‎ (уату... 


نفترض أن хі+1,уі+1‏ فريبة جدا من الحل لذلك 
F1(xi+1, yi+1)=F2(xi+1,yi+1)=0‏ 


E 9‏ كذلك 
Хна Xi ЕВ‏ 
k‏ ۔ Уні‏ 
لحصلنا على الصورة 
=F1(x;y;+ 21+ 51×‏ 0 
ox ду‏ 
2F2(x,y)4- 22+ 92 к‏ 0 
дх ду‏ 
ومنها 


FI ج251 بو‎ - F1G;yp 
0x dy 


2۳2+ of k = - F2(x,y;) 
У 








Ox 
لايساوي‎ Jacobian وهذه المعادلة غير متجانسة ونحصل علی حل عندما محدد‎ 
الصفر‎ 
OF dF 
UP 
J- F% F% 0 
dy dy 
ونحصل على المعاملات كالتالي:‎ 
— Fl(xi, yi) کت‎ 
y 
H= IJ 
— F2(xi, yi) шо 
OF. = ЕК, yi) 
K= /[ 
РА = F2(xi, yi) 
غير‎ х,у نستخدم المعادلتین التالیتین للحصول على قيمتين جدیدتین لكل من‎ 
الافتراضیتین‎ 
0002-32 
Үш=уї+К 


نعيد B SII‏ من جدید باستخدام القيمتين الجدیدتین . 

ثانيا: طريقة نیوتن رافسن المعدلة 

من الواضح صعوبة الطريقة السابقة في حالة تزاید عدد المعادلات (n)‏ أي 
سیکون لدينا (n?)‏ من المشتقات الجزئية لكل معادلة п‏ من المشتقتات 
والمتغیرات , لذلك تعدل الطريقة كالاتي: نحصل على قيمة х‏ الجديدة من طرح 
الحالية من إحدى الدالتین على تفاضلها الجزئي بالنسبة ل × آما قيمة y‏ الجديدة 
فمن طرح الحالية من الدالة الأخرى على تفاضلها الجزتي بالنسبة ل у‏ . 
واختیار أي الدالتین نستخدمه مع × أو y‏ ليس عشوائيا إنما یکون الاختیار 


Lage‏ حیث يؤدي الاختیار الصحیح إلى تقارب نحو Li Jal‏ الخاطین فسوف 


у مع‎ F2 зх مع‎ F1 يبعدك عن الحل. وهنا كمثال اخترنا‎ 
F (xi, yi) 
E (1) 
ox 
. wy. F (xi, yi) 
٢1 OF (ai, yi) ооо фо ооо ооо о ооо ооо ооо ооо ооо фо оо (2) 


дх 


Xj+1=Xj- 











بالنسبة للدالتین : y=x+sin(x)‏ و 3sin(x)-y=0‏ 
واللتان يُراد ЫШ‏ الحل لهما أي نقاط تقاطعهما ويبين الرسم في الصفحة التالية 
نقاط التقاطع. وسنقوم بایجاد Jal‏ بالطریقتین: 





بطريقة نیوتن رافسن 
F1(x,y)=x+sin(x)-y‏ 
F2(x,y)=3sin(x)-y‏ 
نفاضل كلا الدالتين بالنسبة x)‏ وبالنسبة yd‏ 
9F1_1+cos(x) 9F |‏ 
ox oy‏ 
oF? =3cos(x) BLA‏ 
Ox ду‏ 
1+cos(x) -1‏ 
J= = -(1+cos(x))+3cos(x)‏ 
3cosk) —1‏ 
паній (i)‏ ([2008 = 
(x t sin(x) — у) —1‏ — 
H- /J = ((x-rsin(x)-V)-(3sin(x)-9)//J‏ 
(3sin(x)— y) —1‏ — 
kie ika (ii)‏ وا اا و SG‏ 
cos(x) — (x sin(x)— у)‏ +1 
K- IJ-‏ 








3cos(x) — (3sin(x) - у) 


=(-(3sin(x)-y)(1+cos(x))+3cos(x)(x+sin(x)-y))/J 

=(- (x)+y-3sin(x)cos(x)+ycos(x)+3xcos(x)+3sin(x)cos(x)- 
3ycos(x))/J 

=(-3s Se اھ‎ 

= ( y-3sin(x)+cos(x)( 3x-2y ))/J .. — i) 









: X, y 

:Xnew, упем‏ تمثل النقطة الجديدة 

عه : Jules‏ جاكوبي ويساوي المعادلة (i)‏ 

(ii) ويساوي المعادلة‎ хо معامل تفاضل الدالة بالنسبة‎ ІН 

ik‏ معامل تفاضل الدالة بالنسبة لس ویساو d.‏ المعادلة (ددد) 
قدار الخطأ المسموح به وقد اعتبرته 0.00001 












هل jac=0‏ إذا كان نعم.... اطبع ال 

k-(y-3*sin(x)*cos(x)(3*x-2*y))/jac و‎ h=(x-2*sin(x))/jac اجعل‎ 
упеу=у+К و‎ xnew=x+h اجعل‎ 

9. أضف 1 إلى i‏ 

( ا. اطبع ы y e‏ یہ 











1. هل القيمة المطلقة للفرق بين × الحالية والجديدة Jš)‏ من E‏ إذا كان 
لا..... استبدل × الحالية بالجدیدة انتقل إلى 5 
2. هل القيمة المطلقة للفرق بين у‏ الحالية والجديدة Ж‏ من E‏ إذا كان 


لا.... استبدل y‏ الحالية بالجديدة انتقل إلى 5 
13 النهاية 
المخطط الانسيابي 


Flow Chart 
mara ws 
= کر شس و‎ 


X=xnew 


= ° 
Tar-—) -1 
yes 
no : 


H=(x-2*sin(x))/jac 


K=(y-3*sin(x)- 


Divide by zero 
Xnew=x+h Tacohin-ü 


WVnaow- vl 


۶ 
“С 





уе$ 


1 
ves 





f include<conio.h> 
f include<stdio.h> 
f include<math.h> 
void main() 
{ FILE stream; 
int i; 
float x,v,xnew,vnew, jac,h,k; 
printf ("Enter initial value of x, у :"); 
scanf("%f %f", &xnew, &ynew) ; 
i=0; 
stream = fopen("newraf.FIL", "wt"); 
fprintf (stream, "The initial value of x,y is (%2.1f 
, 52۰1۴ ) \n",xnew, ynew) ; 
orint вна, "ба ЗЛЕ ЕЕ УЕ. 
fprintfstteam,"n  . . .— .- . h—  .—  —.— - MES 








fprintf(stream,"\n %d\t %f\t %f\t",i,xnew, ynew) ; 
do 
{ x-xnew; 
y=y ۱ 
јас=2*соѕ (х) -1 ; 
if( јас ==0) { 
fprintf(stream," divide by zero , Jacobin=0 "); 
break; ) 
hz (x-2*sin (x) ) / jac; 
k= 9-070007709779 






۰ 
7 








اسم 
itt;‏ 
fprintf (stream, "\п %d\t $fit $SfNt",i,xnew,ynew);‏ 

}while (fabs (х-хпем)».00001 ع‎ fabs (у-упем)>.00001); 
fclose (stream); 


} 


الادخالات و النتانج 





- 


х,у ай‏ الابتد اتیات 








Enter initial value of x, y :2 1 


The initial value of x,y is (2.0,1.0) 


i x 





ьо Q N P. 


Enter initial value of x, y 
The initial value of x,y 


+ x 


-2 


&UNHOI| 


Enter initial value of x, 
The intial value of x,y is 


i x 


-1. 
-1. 
-1 


1.900996 
1.895512 
1.895494 
1.895494 


.000000 


-1.900996 
895512 
895494 
.895494 





0 
1 
2 
3 
4 


= 0 


-0. 


0.500000 
.107617 
0.000840 
000000 
0.000000 


بر ده نح رمح دح 


0.500000 - 





.000000 
.851493 
.843267 
.843241 
.843241 


الجذر الثاني 
el‏ = 


is (-2.0,-1.0) 





و 





у 20.5 


(0,5,0,5) 





-0.161425 
0.001259 
-0.000000 
0.000000 


آما الطريقة المعدلة 
لو أخذنا الدالتین بهذا الترتیب 


FI (x,y)=3sin(x)-y 
F2 (x,y)=x+sin(x)-y 


ونأخذ للأولى تفاضل بالنسبة ل × والثانية بالنسبة ل y‏ 


arı 


=3cos(x) ed 


ду 


Хил = хі - 3sin(xi)-yi)/ (3cos(xi)) 
Yin = yi — (xi+sin(xi)-yi)/(-1) 


أبدا 

ادخل قيمة х,у‏ الأولية. 

اجعل 1-0 وهي تمثل عدد التكرارات 

اطبع х,у А‏ , 1 الأولية. 

ynew=y-(x+sin(x)-y)/(-1) s xnew=x-(3*sin(x)-y)/(3*cos(x)) اجعل‎ 
الجديدة‎ Lxy ad التي تمثل‎ 1 , xnew, ynew as اطبع‎ 

هل القيمة المطلقة للفرق بين × الحالية والجديدة أقل من E‏ إذا كان 
لا..... استبدل x‏ الحالية بالجديدة (xnew)‏ ثم انتقل إلى 5 

هل القيمة المطلقة للفرق بين у‏ الحالية والجديدة أقل من V Е‏ كان 
لا..... استبدل у‏ الحالية بالجديدة (ynew)‏ ثم انتقل إلى 5 
النهاية 








ho kh ЕБ 


а 


.10 


لمخطط الانسيابي 
flow chart‏ 
| مم | 


X-xnew , 


Xnew=x-(3*sin(x)-y)/(3*cos(x)) 


Ynew=y-(x+sin(x)-y)/(-1) 


Pe 
سرت‎ 





1-1-1 


11 


ves 
no 
ves 


End 
(5-2 (شکل‎ 





#include<conio.h> 

# include<stdio.h> 

# include<math.h> 

void main() 

{ FILE *stream; 

float х, у, xnew, ynew; 

printf('Enter initial value of x, y"); 
scanf ("Sf bf", &xnew, &ynew) ; 
i=0; 

stream = fopen (''modnera.FIL', "м+"); 















пЁ (stream, "The initial value of x,y is 
($2.1f ,%2.1£) \n", xnew, ynew) ; 

fprintf (stream, "\п i\t x\t\t уе"); 
fprintf(stream,''in 





о -‏ سس تس سس 


fprintf (stream, "\п $d ۶ 
$f',i,xnew,vnew); 
{ нні 





(oos 5;‏ ار وود تح سودق سروس ور 

ynew-y- (x+sin (х)-у) / (-1) ; 
itt; 

fprintf(stream,"\n $d $£ 

a 











ae سو‎ ри 
fclose (stream); 


Enter initial value of x, y 





The initial value of x,y is (2.0,1.0) 


л سر ده ين ھ‎ Фф 


14 


18 


20 0 
21 
22 
23 


prone 
1.757074 
1.697342 
2.0 سیق‎ 


' 811609 
.813539 
.992840 
.954219 
.844689 
.852782 
.941547 
.925796 
.867622 
.873738 
.917870 
. 910927 
.880843 
.884470 





NNNNNNNNNNNNNDN N 


.965549 
.912052 
.782754 
„784222 
. 905093 
.881608 
.807414 
.813287 
.873602 
.863443 
.823891 
.828201 
.858242 





'833162 





2. 760836 
2.735107 


43 


45 
46 
47 
48 


49 


50 
51 
52 
53 
54 
55 
56 
57 
58 


59 


1.906570 
1.903323 
1.887941 
1.901036 
1.899461 
1.891635 


1.897503 
1.893532 
1.894064 
1.896901 
1.896511 
1.894499 
1.894770 
1.896206 
1.896009 
7ھ 


1. ' 895 55 
1.895239 
1.895309 
1.895676 
1.895626 
1.895365 
1.895400 
1.895586 
1.895561 
1.895429 
1.895447 
1.895541 
1.895528 
1.895461 


2.835676 
2. 850726 
2.848544 
.838070 
.839426 
.847000 
845935 
840606 
.841314 
.845137 
.844607 
.841903 
.842267 
.844199 
843934 
842563 
.842748 
.843725 
.843592 
.842898 
.842992 
843486 
.843067 
.843115 
.843365 
.843331 
.843153 
.843177 
.843304 
843287 
2.843197 
2.843209 
2.843273 
2.843264 


5 


05 


е 


° 


. 


е 


. 


Enter initial value of x, 


у 


:0.5 


0. 


5 


The initial value of x,v is (0.5,0.5) 


i 


о بن لطر كن وہ لك‎ N شر‎ O 


ooooooooo, 


سای .0 





0. ‚217908 
0.056940 
0.144872 
0.037329 
.074650 
.024699 
.064273 
.016411 
042838 





` 019036 
00 





'003233 
.008460 
.002155 
.005640 
0.001437 


6 ۵ © © © © 


سیت 





“001671 
. 000426 


о 


у 


0.500000 
0.979426 
0.286733 
0.651855 
0.177078 
0.434095 


0.113849 


0.289237 
0.264690 
0.192775 
0.049396 
0.128502 
0.032821 
0.085663 
0.021848 
0.057107 
0.014556 
0.009701 
0.025381 
0.006466 
0.016920 
0.004311 


0.002874 
0.007520 
0.001916 
0.005013 
0.001277 
0.003342 


.001114 
.000284 
.000743 
.000189 
.000495 
.000126 


.000851 
.002228 
. 000568 
.001485 
.000378 
. 000990 
.000252 
.000660 
. 000168 
.000440 
.000112 
.000293 
. 000075 
.000196 
. 000050 
.000130 
. 000033 
. 000087 
. 000022 
. 000058 
. 000015 
. 000039 
. 000010 





. 000084 
.000220 
. 000056 
. 000147 
.000037 
.000098 
.000025 
.000065 
. 000017 
. 000043 
: РО 


6 © © © 





A ۱ 
о 


© © © © 


`000007 
. 000019 


0.000005 
_. 000013 


ло 
о 


(л 
N 
о 


۳ 


nter initial value of x, у :-2 -1 


The initial value of x,y is (-2.0,-1.0) 





3. | -2.909297 
- 137745 -3.143961 
024 -2 9 
.697342 a 4 
-2.321275 -2.689346 


1.000000- مہہ ہو 








Q در نم ين صر‎ о 
۱ 
= 


57 -1.895541 -2.843209 
58 -1.895528 -2.843273 


59  -1.895461 


-2.843264 






у V الجديدة أي‎ 


ре من‎ 
Хы = xi - (3sin(x))-vi )/ (3*cos(xi)) 
Yin = yi - (Xn FSIN(X,,1)-VI)/(-1) 





З *вїп (х ) -y) / (3*cos (x)) ;‏ يوه 
ynew-y- (xnewrsin (хпем) -у) / (-1) ;‏ 


فکانت النتانج لنفس الادخالات کالتالي: 
Enter initial value of x,v 1‏ 





The initial value of x,v is (2.0,1.0) 


2.000000 


. 000000 


3.384041 
2.057167 
1.851002 
1.872380 
1.914300 
1.883480 
1.904657 
1.889309 
1.900055 
1.892333 
1.897785 
1.893886 


143961. 
941202. 
812001. 
870419. 
۶ھ 
834992. 
849441. 
2.839011 
2.846337 
2.841084 
2.844799 
2.842145 


2.844028 1.896649 14 
2.842685 1.894678 15 
2.843639 1.896078 16 
2.842959 1.895080 17 
2.843442 1.895790 18 
2.843099 1.895285 19 
2.843343 1.895644 20 
2.843169 1.895388 21 
2.843293 1.895570 22 
2.843205 1.895441 23 
2.843267 1.895532 24 
2.843223 1.895467 25 
2.843255 1.895514 26 
2.843232 1.895481 27 
2.843248 1.895504 28 
2.843237 1.895487 29 
2.843245 1.895499 30 
(جدول 3-2( 


Enter initial value of x, y :0.5 0.5 


The initial value of x,y is (0.5,0.5) 


I x y 

0 0.500000 0.500000 
1 0.143613 0.286733 
2 0.095576 0.191007 
3 0.063669 0.127295 
4 0.042432 0.084850 
5 0.028283 0.056563 
6 0.018854 0.037708 
7 0.012569 0.025138 
8 0.008379 0.016759 
9 0.005586 0.011172 
10 0.003724 0.007448 
11 0.002483 0.004966 
12 0.001655 0.003310 
13 0.001103 0.002207 


.001471 
.000981 
.000654 
.000436 
. 000291 
. 000194 
. 000129 
. 000086 
. 000057 
. 000038 

(جدول 2 -4( 

و بتغییر اختیار الدوال(إعادة الترتیب) كالاتي: 
F1(x,y)= x+sin(x)-y : F2(x,y)= 3sin(x)-y‏ 

ونأخذ للأولى تفاضل بالنسبة ل x‏ والثانية بالنسبة ل y‏ 

ے77 
ду‏ 


14 
15 
16 
17 
18 
19 
20 
21 
22 
23 


.000736 
.000490 
.000327 
. 000218 
. 000145 
. 000097 
. 000065 
. 000043 
. 000029 
. 000019 


Oo‏ ہ ہ ہ ہ ہ ہ ه ه 
ه مه مه ہ ہ و ہہ 


arı 


=1+cos(x) :‏ 
ox‏ 
ونطبق القانون 
Хил = Xi - (xi+sin(xi)- yi)/ (1+cos(xi))‏ 
Үш = yi - Gsin(xi)-yi)/(-1)‏ 
فان الناتج سیکون مبتعدا عن الحل وهذا هو ما حصل فعلا 


0 3.000000 2.000000 
1 -111.026443 0.423360 





71 -1431034113556480.000000 -1.920399 
72 49127493789024256.000000 -0.708711 
73 +NAN +NAN 


الاستنتاج 

* في الطريقة المعدلة وبعد استخدام قيمة x‏ الجديدة у slay‏ الجديدة تقلص 
sas‏ التکرار إلى النصف تقریبا عن الطريقة المعدلة التي استخدمت فیها قيمة x‏ 
القديمة لایجاد у‏ الجديدة أي وجدنا الحل في نصف الوقت, 

* في طريقة نیوتن رافسون تم الوصول إلى الحل بأقل تکرار من المعدلة 
وبتقارب دائم نحو الحل Lain‏ الطريقة المعدلة نجدها مرة تتقارب وأخرى BIS‏ 
مما يزيد عدد التکرار. وهذا یتضح من الأشكال التوضيحية لاحقا. 

Ж‏ اختیار أي معادلة نعتبرها FI (x,y)‏ مهم. حيث نلاحظ التقارب والتباعد الناتج 
عن ذلك 

* نلاحظ من الجدول, اختلاف نتائج الطرق سالفة الذكر من حيث عدد التكرارات 
حسب قيمة х,у‏ المدخلة مع ملاحظة أي من الحلول الثلاثة نصل إليه في كل 
مرة. ومن المعلوم أن هذه التكرارات تعتمد على E‏ (قيمة الخطأ المسموح به) . 






































نيوتن رافسن المعدلة مع | نيوتن رافسن . . . . 
l .‏ " نیوتن رافسن | 5 
استغلال x‏ الجديدة لایجاد y‏ المعدله النقطة المدخلة 
التکرار < الحل التکرار > الحل | التکرار > Jal‏ 
з ۳ (0,0) > 4‏ 2.8)>5-,1.8-( | 5-,3- 
4 > )0,0( 4 > )0,0( 5 < )0,0( 0.6,0.6 
.1.8( 
2.8)>33-,1.8-( 2.8)>75 1.8,2.8)>5( | 1.5,0.75 
0,0)>27( تباعد 1.8,2.8)>5( | 1.5,2.5 
1.8,2.8)>29( 1.8,2.8)>57( | 1.8,2.8)>4( 2,2 
0,0)>24( 1.8,2.8)>61( | 1.8,2.8)>5( 3,5 
s% (0,0)>46‏ )1.8,2.8)>114(| 74 
0,0)>30( تباعد 2.8)>7-,1.8-( 10,12 
0,0)>31( تباعد 1.8,2.8)>145( 20,20 














(جدول 2 -5) 





شكل توضيحي یبین مقارنة بین الطرق من حیث التکرارات 


مقارنة بين عدد التکرارات لایجاد الحل بطرق مختلفة 


10 


60 


50 





40 


30 


20 











2 ن‌رم‎ laso 





(شکل 2 -6( 
ن.ر : طريقة نیوتن رافسون 
ن.ر.م 1 : طريقة نیوتن رافسون Ала)‏ 
ن.ر.م 2 طريقة نیوتن رافسون المعدلة والتي قمنا فیها باستغلال قيم × الجديدة 
لایجاد y‏ 





شكل یوضح نقارب 
ال Мыз‏ )2,1( إلى 
Jal‏ في طريقة 
نیوتن ر افسن 











شكل يوضح تقارب 
النقطة (2,1) إلى الحل 
في طريقة نيوتن رافسن 
المعدلة 





(شكل 2 -8) 








3 
2 شكل یوضح تقارب 
s-1(2,1) ТЕ,‏ 
1 الحل في طريقة 
Р‏ نيوتن رافسن المعدلة 


0 1 2 3 4 الجديدة لایجاد y‏ 











(شكل 9-2( 
الأشكال الثلاثة في شکل واحد مع بیان نقطة البدء ونقطة الحل 





خط یوضح نقارب 
النقطة )2.1( إلى 
JA‏ في طريقة 
نيوتن رافسن المعدلة 
الجديدة لإيجاد у‏ 





N 05 0‏ 15 2 25 3 35 4 النقطة (2,1) إلى 
ال i ы‏ 











النقطة )2,1( إلى 
الحل في طريقة 
نيوتن رافسن 


4-2 :43 الموقع الخاطی 
False position‏ 
зь Да в TS E I‏ 
نظرا لأهمية هذه الطرق ¢ رأيت أن أضيف هذه الطریقة و التي لا تقل آهمية عن 
الطرق السابقة » وللمزيد من الاتساع في الموضوع و لان التعدد يفتح آفاق الخيار 
امام الدارس من اجل التنوع و الاستفادة الشاملةء كنت اريك أن ابدا في شرح ) 


Y (Lagrange interpolation polynomial‏ اعرف ترجمة Аз ус‏ دقيقة لها؛ 
على کل حال لنبدا في طريقة ( الموقع الخاطی) ) (false position‏ 


المفهوم النظري 
إذا كانت الدالة معرفة في الفترة ) x2‏ , 61 و كان f(x1)*f(x2)«0‏ 
من الرسم الموضع یمکن استنتاج أن المثلئین : 

О, f(x1) ,f(x2) 

×3, ×2 , (X2) 
: متشابهان فیکون‎ 
F(x2)/(f(x2) — f(x1)}=(x2 — x1) (x2 — x1) 
×3 و تکون قيمة‎ 
X3=x2 — f(x2)* ((x2- x1)/f(x2)-f(x2) } 
و يمكن كتابة الخوارزمية للمعادلة التالية‎ 

у1=х**3 +4*x**2-10 

مع ملاحظة أن الرمز ** يعني في لغة الفورتران للقوة الثانية مثلا أو الثالثة أو 
غيرها أي الأس 
ors‏ 


start .1 © 
۲6۵2001 ,х2) .2 e 


do .3 e 
yl=x1**3 +4*x1**2-10 6 
y2=x2**3+4*x2**2-10 e 
X3=x2 — f(x2)*((x2- x1)/f(x2)-f(x2)) ° 
Y3=x3**2+4*x3**2-10 ° 
if(y3<0.000001)goto (8) 
if(y1*y320.0)goto (7) 
x]=x1 6 e 
X3=x2 ° 
goto (3) ° 
хівхі .7 ° 
X232 ° 
goto 3 ° 
write x3,y3 .8 e 
stop .9 e 


و البرنامج بلغة فورتران كالاتي : 


Read(*,*)x1 x2 e 
ylzx1**3 +4*x1**2-10 e 
y2=x2**3+4*x2**2-10 e 
X3=x2 — f(x2)*((x2- x1)/f(x2)-f(x2)) ° 
Y3=x3**2+4*x3**2-]0 e 
if(abs(y3).1t.0.000001)goto Š ° 
if(y1*y3.gt.0.0)goto 
xlexl 
x3=x2 
goto 3 
xisxi 
x2=x2 
goto 3 
write(*,*) x3 ,y3 
stop 
End 


л B 
e e 





(شکل 11-2( 
و بهذا نكون قد انهينا طريقة الموقع الخاطی 


الباب الثالث 


الحلول العددية للمعادلات التفاضلية العادية 


3 -1 الحلول العددية للمعادلات التفاضلية العادية 


سوف نقوم بحل المعادلات من النوع 
dy/dx = f(x,y)‏ 


أي معادلة خطية من الرتبة الأولى . ولحل هذا sill‏ £ من المعادلات هناك 
مجموعة من الطرق منها: 

Euler's Method 19 طريقة‎ 2> З 

نفترض أن الحل y=F(x)‏ متصل وقابل للتفاضل وأن نقطة البداية هي 
(x0,y0)‏ وأن الحل مطلوبا عند ×=« لذلك يجب تجزئة Дай‏ إلى 0 من 
الاجزاء برض ule ( w‏ أن تکون صغیرة) 


w =(x-x)/n cs 


ثم نطبق القانون 
yii=yi +w * f(xi,vi)‏ 
I = 0,1,..n-1 Gus‏ 
ملاحظة: 
- یکون الخطأ آقل بقيم WI Ji‏ 
- تزداد الحسایات ola jb‏ عدد да‏ ات آي تصغیر W‏ 
US -‏ ل تعتمد على التي قبلها مما بنشر الخطأ في حالة وقوعه. 
امتداد طريقة 9№ Extended Euler Method‏ 
بأخذ ثلائة حدود من متسلسلة تایلور التي حصلنا منها على EM‏ نحصل على 
EEM‏ 
Урагуї +w * f(xi,yi) + w2/21 * f '(xi,vi)‏ 
I-0,1,.n-1 Cus‏ 
3-3 طريقة اویلر الاكثر امتدادا More Extended Euler‏ 
Method‏ 
بأخذ الحد الرابع من متسلسلة تایلور dax‏ على MEEM‏ 
yii=yi +w З f(xi,yi) + w2/21 * f `(xi,yi) + ми /31 * f ``(xi,yi)‏ 


Modified Euler Method طريقة اویلر المعدلة‎ 4- 3 


في هذه الطريقة نستغل EM‏ لایجاد 1+1 7 (Predictor)‏ ثم نستخدمها 
аа‏ 
avg- ( f(xi,yi) + #(хі+1,у`і+1) 2‏ 
نحسب من (Corrector) уі+1 iss‏ 
у1+1 = yi + w* avg‏ 


5-3 طريقة رنج کوتا Runge-Kutta Method‏ 


تتميز بدقتها وجودنها 
Yi+1= vi + (k1+2k2+2k3+k4)/6‏ 


КІ = w* f(xi,yi) 

K2- w* f(xi+w/2,yi+k1/2) 
КЗ = w* f(xi+w/2,yit+k2/2) 
КА = w* f(xi+w,yit+k3) 


وبأخذ السؤال الذي يقول آوجد حل مسألة القیم الذاتية ؛ 
dy/dx =—ху ; x0=0 , Wel‏ 
وذلك عند 1= ¿Ü Lle x‏ ۷-0.1 (قارن بقيم (2/×-)۸) 
سنقوم بحله بكل الطرق المذكورة آنفا 
f(x,y)=-xy Cus‏ 
اویلر 
Yi+1 = yi - w * xi*yi‏ 
بتفاضل الدالة Ë‏ وزيادة as‏ للسابقة نحصل على امتداد اويل 
У`=-ху f y`=-xy`+t(-1) = -ху зу‏ 
بالتعويض عن قيمة y‏ 
Y" = -x(-xy)-y = x у-ує У(х -1)‏ 
allis,‏ نحصل على امنداد اویلر 
Yi+1 = yi - w * xi*yi + w^*yi*(xi^-1)/2‏ 
بتفاضل f‏ مرة آخری وزيادة حد نحصل على اویلر الاکثر امتداد 
y'-yxX-Dt у``=у(2х)+(х^-1)у`‏ 





y Ax بالتعویض عن‎ 


على اویلر الاکثر امتدادا 








Y 'I+1=yi —w*xi*yi 
Avg= ( -xi* yi - xi+1*y i+1)/2 
Yi+1 = yi + w* avg 








Virl- yi + ( (Е1+2К242К3+К4)у6 


يقة نكرر العمل عدد n‏ من المرات بعدد الفترات المأخوذة . 








АЧ x عند‎ y قيمة‎ ۷ 


:K1, k2, КЗ, 4‏ معاملات رنج کو 








#: عرض افترة الوا 
Wall Аа ۴‏ ( التفاضل) 












6 ادخل عدد الفترات n‏ وقيمة х,у‏ الابتدائیات 
8. احسب F=-x*y‏ 

kl=w*F сыз! 19 

F=-(x+w/2)*(y+k1/2) Gus! .20 
F=-(x+w/2)*(y+k2/2) احسب‎ .2 

k3=w* F احسب‎ .23 

F=-(x+w)*(y+k3) احسب‎ .24 

k4-w*F сыз!) .25 
ynew=y+(k1+2*k2+2*k3+k4)/6 احسب‎ .26 








اجعل y=ynew ,x-xnew‏ 
9. اضف 1 إلى i‏ (121+1) 
I=n 08.30‏ إذا كان لا ارجع إلى 5 
1. اطبع ynew Дай‏ 






المخطط الانسيابي 
Flow Chart‏ 





ا هم | 


F=-(x+w/2)*(y+k1/2) 
120-7 








F=-(x+w/2)*(y+k2/2) 
K3=w*F 


F=-(x+w)*(y+k3) 
K4=w*F 


Ynew=y+(k1+2*k2+2*k3+k4)/6 
Xnew=x+w 


End 
(1-3 (شکل‎ 
‚ ملاحظة‎ 

- البرنامج الاول هو البرنامج الرئيسي الذي یستخدم طريقة رنج US‏ 
بالعرض المطلوب في السؤال w=0.1‏ 
- البرنامج الثاني جعلت فيه W‏ متغيرة بحيث تم أخذ فترات مختلفة 
- البرنامج الثالث استخدمت طريقة اويلر وإمتداداتها 
- البرنامج الرابع طريقة اويلر المعدلة 
- من الممكن طبعا ادماج الطرق جميعا في برنامج واحد 


التعلييق 
- طريقة رنج كوتا هي أدق الطرق 
- اويلر المعدلة اعطت نتائج ادق من اويلر الاخريات 
- أما اويلر وامتداداتها فان الاكثر امتدادا كانت أدق ثم امتداد اويلر وأخيرا 
اويلر 


- بزيادة عدد التكرار أي بتصغير w‏ تكون النتائج أدق 


06055 > 





(شکل 2-3( 
العلاقة بين مقدار у Аай у уу‏ الناتجة Cus‏ یتضح التقارب نحو الحل الصحیح 
بتقلیل да)‏ ات 





(شکل 3-3( 


المحور السيني یمثل قيمة X‏ من 0 إلى 1 بزيادة قيمة W‏ والمحور الصادي هو قيمة 
y‏ 
هذا الشکل ينيم الطرق ناک 





(شکل 4-3( 


منحنی الدالة الأصلي مع القیم الناتجة من طريقة رنج US‏ واویلر المعدلة 


6-3 طريقة الرمي 


مقدمه 

تستخدم طريقة الرمي لایجاد حل المعادلات التفاضلية من الرتبة الثانية , ولذلك 

يجب علينا معرفة ثابتين وعادة ما تکون الثوابت Le}‏ الدالة ومشتقتها الأولى 

عند نقطة البداية (مسألة القيمة الابتدائیة) أو الدالة ومشتقتها عند نقطتين 

مختلفتين (مسألة القيم الحدية) كمسألتنا هذه. 

وتتلخص طريقة الرمي في: 

- نفترض الشرط الناقص من شروط مسائل القيم الابتدائية وهو هنا تفاضل 

الدالة عند x=]‏ 

اوجد الحل بالافتراضات الجديدة )73 (x‏ وقارنه مع الشرط المعطى عند 

- آعد تغییر القيم الابتدائية أي الفرض الذي فرضته حتى نحصل على 
المطلوب وهو (في هذه المسألة y(3)=10.0179‏ ( 


وبأخذ السؤال الذي يقول أوجد حل مسألة القيم الذاتية : 
=y ; y(1) =1.1752 , y(3) 210.0179‏ 06/0۰ 
آولا: نأخذ قيمة افتراضية ل ۷ عند 1=× ثم نستخدم طريقة رنج US‏ لایجاد قيمة у‏ 
х=3 AM‏ ) 
ثانیا : نفترض قيمة آخری ل ۷ عند 8-1 ونعید استخدام طريقة رنج $ U‏ انحصل 
على y‏ عند «же?‏ 
ثالثا نستخدم التوليد الخطي لإيجاد قيمة y^‏ أخرى كالتالي: 
Y X VY‏ 
RI gl 22-1‏ 
R2 g2‏ 
D xo[2]‏ 


Y=y0+pV yO 
P=(xp-x0)/h 


با el o‏ مھ T‏ اتنا 


H-r2-ri 

717-02-1 

P=(d-r1)/(r2-r1) 

Y= а вы (d-r1)/(r2-r1)*(g2-g1) 





x das: tstart‏ ( الابتدائیة) 





:xstart‏ قيمة у‏ الابتدائية 





м بمقدار‎ ух 448 متغیر يأخذ‎ БО 

















كنا قد تناولنا في المقالات السابقة» طرق حل : معادلة ~= ذات مجهول 
sala‏ فقط» وسر أن 0 من if) ы"‏ | الثانية فماذ 





l-4‏ طريقة جاكوبي لحل مسائل القیم الذاتية 


مسائل القیم الذاتية تتلخص في ایجاد À ad‏ وهي القیم الذاتية و X‏ وهي المتجهات 
АРЕ‏ رة зів‏ 

Ax= ХІХ 

all al2 3 хі м 0 0 хі 


a22 а23 |*| х2 |=| 0 А 0 |*| x2‏ 21و 
a31 a32 3 х3 0 0 23 х3‏ 


حيث A‏ مصفوفة مربعة متجانسة 
طریقه جاكوبي : 
بهذه الطريقة ay ys‏ الوصول إلى المعادل 


Bx= AX 
مصفوفتان فطریتان‎ В و‎ А Cus 
bll 0 0 
B= 0 b22 0 
0 О 33ط‎ 
ii < A JENS ذات‎ да . e PAR jiġ 
وطريقة جاكوبي تقوم بدوران المصفوفة بزاوية 6 في مستوی آخر‎ 
Хх 


T cus‏ مصفوفة الدوران ولأنها في ثلاثة ابعاد توجد منها ثلائة آنواع 


cos0  -sinO 0 
Tl= 500 وم‎ 0 
0 0 1 
cos0 - $10 
T2= 0 1 0 
ѕіп Ө 0 соѕӨ 
1 0 0 
T3= 0 056  -sinO 
0 sin 060 








بما آن 


AX-AX 
Atx-ATX 
T أن ضرب المصفوفة في معکوسها =1 لذلك نضرب الطرفین في معکوس‎ Lo 
T'ATx-AT ‘Тх=Ах 
نختار الزاوية المناسبة أي أن يكون الحد‎ B مصفوفة فطرية‎ ТАТ ولكي نجعل‎ 
الغیر الطري في المفوفة الناتجة من الضرب يساوي الصفر‎ 
al2í(cos2 O-sin2 Ө) +соБ Osin 9(a22-a11) 
tan (29) =2*а12/ (а11-а22) 
Ө=+ап-1 (2*а12/ (а11-а22) ) /2......(1) 
وللحصول على المتجهات الذاتية نقوم بضرب مصفوفات الدورات المستخدمة من‎ 
AM البداية وحتی الحصول على القيم‎ 
Verl*T2*...TIH 


ca E الدو‎ Adm pp d ца 


في المثال المعني 

2 -1 0 
-1 2 -1 
0 -1 2 


Àz 








1. نأخذ أكبر عنصر غير قطري في المصفوفة وهو هنا 1 - 

2. مما سبق نأخذ TI LJ‏ أو ОУ ТЗ‏ أكبر „ае‏ مکرر في مکانین 

3. نحسب فيمة الزاوية من المعادلة (1) وتساوي 

O-tan-1(-2/(2-2)—-1/2 

4. نعین المصفوفة الناتجة من الضرب T AT‏ وتصبح هي А‏ 
نضرب مصفوفة الوحدة في T‏ المرة الأولى ثم الناتج في T‏ المستخدمة 
أخيرا. وذلك للحصول على المتجهات AGA‏ 

6 نعید الخطوات 1و2و3 على المصفوفة الجديدة حتی نصل إلى أن تكون 

















:BigerO‏ تحدد آکبر عن 
مت | تساوي ا 

















Tetha—PT/ Tetha—PT/ 


ttii=ttjj=tii=tjj=cos(te 
tha) 


A=tt* A*t 


al2&al 


29.21- 





#include<stdio.h> 
#include<conio.h> 
#include<stdlib.h> 


#include<math.h> 

#include<iostream.h> 

#define max 3 

Ant 1.3; ,k,n; 

void equal (float s[][3],float c[][3]); 
void input mat (float a[max] [max]) ; 
void mull_mat (float aa[] [max], float 
b[] [max], float c[] [max]); 

void write mat (float а[] [max]) ; 

void biger (float а[] [max]); 

void tandinvers (float t[] [3], float tt[][3],float 
a[] [max]) ; 

FILE *in,*out; 

main() 

( int k; 

float 

а [тах] [max] , a [max] [max], c [max] [max], tt [3] [3], t [3] 
[3], v[nax] [max]={1,0,0,0,1,0,0,0,1}; 

clrscr(); 
n=3; k=0; 








jacread. fil", "r"); 






ist. 
fclose (іп); 


LP‏ لس وہہ pa E‏ د 










fprint hon 
while( (fabs(a[0][1])>0.0001) || 
(fabs (a[0][2])>0.0001) || (fabs (а[1] [2])>0.0001) 





{ k=k+1; 
biger(a ;( // find biger element of A 
fprintf(out,"\nCycle по. за & The biger element 
t %5.4f\n",k,a[i][j]); 

tandinvers (t, tta 34 / / find Angle , T 
& invers T 

fprintf(out,"T = Win"); 

write mat (t); 

mull mat (tt,a,d); // multiplicat invers T in 












mull mat (d,t,a); // multiplicat А in T 
fprintf(out,"A = Nin"); 

write mat (a); 

mull mat (v,t,c); //clculat Eigenvectors 


equal (у, с); 
) 
ره‎ "The lamda is "); 
[1=0;1<3;1++) { 
E (out, "AnJ$d-$5.4fNn v= ",i,a[il[il); 
for(jz0; j«3; j++) 
fprintf (out, &5,4F ",у[3][1]); 
} 
fclose (out); 
return (0) : 
} 
// function to find angle , T and invers T 
void tandinvers (float t[][3],float tt[] [3], float 
a[] [max] ) 
{ float x,z,y,tunit [3] [3]={1,0,0,0,1,0,0,0,1}; 
if (ali] [i]==a[5] [j1) B x-3.1416/4; 
else x=-3.1416/4; 
else x=(atan(2*a[i] [j]/(ali] [1]-а[51[51)))/2; 
y=cos (x); z=sin (х); 
equal (t, tunit) ; // t= {1,0,0,0,1,0,0,0,1} 
equal(tt,tunit); 
// select TI , T2 or T3 consider biger element 
Y 
for (i=0; i<n; i++) 
+ :3++( 
fscanf (in, "SE", salillji); 
} 
void biger (float a[] [max] ) 
{ 









if (fabs (a[0] [1]) >=fabs (а[0] ]2 ) (1 
if (fabs (a[0][1])>=fabs(a[1][2])) (i-0;j:1;) 
else (i-1;j-2;) 
else if (fabs (а[0] [2])>=fabs(a[1] [21)) 
{i=0; j=2; } 
else {i=1; j=2; } 
1 
void mull mat (float aa[][max],float 
b[] [max], float c[] [max]) 
{int k; 
for (i=0; i<n; i++) 
for (j=0; j<n; j++) 
{ c[i][j]=0; 
for (k=0; k<n; k++) 
c[i][j]=c[i] [3] +аа [1] [k]*b[k] [j] ; 








H 


) 
) 


void write mat(float c[][max]) 


for (i=0; i<n; i++) 


{ 


for (320; j<n; j++) 
fprintf (out, "۲55۰4۶۱۳۲, cli] 7+ 


The Array of A= 


2.0000 
1.0000- 


Cycle по. 1 & The biger element : 


т = 
07072 
04 TOTIS 


Cycle по. 2 & The biger element : 


fa 
0.8881 
0.0000 
0.4597 
А = 

3.9660 


0.0000- 


Cycle no. 3 & The biger element : 


T = 
1.0000 
0.0000 
0.0000 
А = 


1.0000- 
2.0000 


D, TOT 
ول‎ 7071 


0.0000 


б»Фезі- 


D.6280- 


0.0000 
0 +7 
0.5241 


0.0000 
1.0C00= 


0.7071 
2.0000 


0.4597- 
0.0000 


0.8881 


0.0000- 
0 6280- 
1.6340 


0.0000 
0 52416 
0.8517 





-1.0000 


0.7071 


-0.6280 


8ع =0.2768 3.3650 
08 0.6136 -0.2768 
2.0204 0.0000 0.1704 


Cycle no. 4 ع‎ The biger element : -0.2768 
T = 

0,9951 0.0991 0,0000 

0.0991- 0.9951 0.0000 
0.0000 0.0000 1.0000 


44b  б.0000- ء8٤67‎ 
020000 0,5860 00,0109 
0,1695 0.0169 2.0204 


Cvcle no. 5 & The biger element : 0.1695 


0.9927 0.0000 0.1207- 
0.0000 1.0000 0.0000 
0.1207 0.0000 0.9927 
A= 

3.4142 0.0020 0.0000 
0.0020 0.5860 0.0168 
0.0000 0.0168 1.9998 


Сус1е no. 6 & The biger element : 0.0168 
TE 

1.0000 0.0000 0.0000 

0.0000. 0.909 0.0119 

0:8000 0.0119- [9 

A — 

3.4142 0.0020 0.0000 

0.0020 0.5858 0,0000- 


Cycle no. 7 & The biger element : 0.0020 
т = 
1,0000  0.000T7T- 0.0000 


0,0007 1,0000 0,0000 


0.0000 0.0000 1.0000 


3.4142 0.0000 0.0000 


The lamda is 


42 41 . 023 ل 
v= 0.5000 -0.7071 0.5000‏ 


J1=0.5858 
v= 0.5000 0.7071 0.5000 
02-0 
v= -0.7071 -0.0000 0.7071 
الاستنتاج و التعلیق‎ 
وطباعة المصفو 44 ولحساب الضرب‎ bel äl تم استخدام عدة بر امج فرعية‎ ۶ 
وغيره‎ 


- يختلف عدد الدورانات باختلاف الدقة المختارة للخطأ المسموح به فكلما 
زادت الدقة كلما زاد عدد الدوران فالعلاقة طردية (عند 0.0001 عدد 7 
دورانات) 

- نلاحظ في حالة تساوي قيمتين كأكبر قيمة فإن الاختلاف في أخذ T‏ يؤدي 
إلى نفس الحل ولكن ليس بنفس الترتيب . 

- استخدمت شرطا خاصا لإيجاد قيمة الزاوية عندما يكون المقام بصفر 

- لاحظ أن من موقع jb‏ عنصر نستطيع تحدید Т‏ فلو كان موقع أكبر 
عنصر all=cos, akk=cos , alk=-sin , akl-sin а bildl-Ij-k‏ 

- استخدمت مصفوفة الوحدة لوضعها في مصفوفة الدوران قبل تحمیلها كما 


- 


سبق . 


Lagrange's Interpolating Polvnomial جر انج للتو ليد‎ У 44, p^ 2- 4 


ضا طريقة لاجرانج у а‏ جیدة لتولید pd‏ جدیده ارت 
الحدود وک دك لایجاد صيغة JIS‏ 5 الحدود . وهي تصلح لأي عدد من 
النقاط . وتعتمد درجة كثيرة الحدود الناتجة على aae‏ النقاط المستخدمة. 
درجة كثيرة الحدود = 11-1 حيث n‏ عدد النقاط 
فلو كانت عدد النقاط 3 فستكون كثيرة الحدود من الدرجة الثانية وهكذا 
وتكون axe‏ معاملات لاجرانج L‏ تعتمد على عدد النقاط وتكون قيمة L‏ 
على الصورة: 


Li- (х= xi)(x — x2)......(x — xi - i)(x — xi + 1)....... (x — Xn) I 


С) а) а РИ (Xi Xn) nnee 


ححاصل ضرب حاصل طرح قيمة من کل النقاط عدا النقطة الحالية 
حاصل ضرب حاصل طرح النقطة الحالية من کل النقاط عدا نفسها 
ثم لایجاد المعادلة النهائية : 





Y=L{x) = Yus yi 

ومن الملاحظ أنه لا یمکن بسهولة أو قد AD‏ حلها یدویا خاصة عندما 
تکثر axe‏ النقاط ویزداد صعوبة عندما نکون ad‏ النقاط بالأرقام العشرية. 
لذلك تحل باستخدام الحاسوب, وسنبرمج طريقة لاجر انج باستخدام لغة C‏ 


قائمة بالمتغیرات المستخدمة في البرنامج 

usi 7,4‏ اکر زر 

М‏ عدد النقاط المدخلة 

laste у النقطة المراد ایجاد قيمة‎ : Xp 

Xi, yi‏ : النقاط المدخلة 

©: حاصل ضرب طرح قيمة xp‏ من النقاط أي البسط في القانون (I)‏ 








ادخل قيمة хр‏ الس اد ایجاد y‏ عندها 
اجعل قيمة С‏ 


احسب می 3 Е.‏ 


هل انتهت النقاط Men‏ 


رمیا 
















المخطط الانسيابي 
Flow Chart‏ 










ves 
no 


D=D*(xp-xi) 
C-C*(vi_vi\ 


yes 


L=C/D 
Qnm-&€nmaT Ж 





f include<stdio. 


# include<conio.h> 
void main() 
(int i,j,n; 
float L[9],sum,xp,x[9],y[9],C,D; 
sum=0; 
printf("Enter n value\n "); 
scanf ("%d", &n) ; 
printf("\n Enter x and у Mn "); 
for (i=0; i<n; i++) 
scanf ("Sf S£",&x[i],&y[i]); 
printf("\n Enter point of x i E 
scanf("$f",&xp); 
clrscr(); 
printf("Lagrange method\n\n") ; 
printf ("\n x ")) 
for(i=0;i<n;i++) Printf("%4.2fNt",x[i]); 
printf("Nn y "); 
for(i=0;i<n;i++) Printf("%4.2fNt",y[i]); 
printf(" хр= %4.2f\n",xp); 
BIIN سے چم سے‎ o مسجم‎ ооо 
for (i=0;i<n; i++) 
{ C=1; D=1; 
for (j=0; ј<п; j++) 
{ if(il-3) ( 
C=C* (xp-x[j]) ; 
D=D* (x [i]-x[j]); 
L[i]=C/D; 
sum=sum+L [i] *y[i]; 
printf ("\п L(%d)= $f ",i,L[i]); 
} 
printf ("Хом 
printf ("\nL=L1*y1+12*y2+....\n"); 
printf 
n\n"); 
printf ("L= $f ",sum); 




















1 Jù 
Lagrange method 








x 0.00 1.00 2.00 3.00 4.00 
у 0.00 2.00 8.00 18.00 32.00 xp= 
5 - 0 0 


L(0)= 1.000000 


L(1)= -5.000000 


L(2)= 10.000000 


L-Ll*yl4l2*y24.... 





L- 50.000000 


تکون النتائج 








x 0,00 30.00 60.00 930,00 


у 0.00 0.50 0.8660254038 1.00 хр= 


L(0)= -0.061728 
L(1)= 0.740741 
L(2)= 0.370370 
L(3)= -0.049383 


L=L1*y1+12*y2+.... 


L= 0.641738 





INN IN NN 000 


Lagrange method 





x 90.00 30.00 0.00 60.00 
у 1.00 0.99 ИНО 0.8660254038 хр= 
40.00 


L(0)= -0.049383 
L(1)= 0.740741 
L(2)= -0.061728 
L(3)= 0.370370 


L=L1*yl+12*y2+.... 


L= 0.641738 





Lagrange method 


x 0.00 30.00 60.00 


40.00 


L(0)= -0.111111 
L(1)= 0.888889 


L(2)= 0.222222 


L=L1*y1+12*y2+.... 


L= 0.636895 


is بنقاط أخرى ولکن‎ 4 2 JE 





Lagrange method 


30.00 60.00 90.00 
0.50 0.8660254038 1.00 xp= 40.00 


L(0)= 0.555556 
L(1)= 0.555556 
L(2)= -0.111111 


L=L1*y1+12*y2+.... 


L= 0.647792 





$ .ففي 2- »© كانت 0× قريبة من النقطة الثالثة‎ a 





كانت xp‏ قريبة من النقطة الأولى في 4-2 زادت L‏ فهي 





Dimension x(50),f950) 
Write (*,*)'enter value of polynomial’ 
Read (*,*) n 
Write (*,*)'enter point of x that you find solution’ 
Read(*,*)x0 
Do 10,i=1,n+1 
Write (*,*)'x (i), «fC, =" 
Кеаа(*,*)х(1), 1) 
10 Continue 
Xp=0 
Do 20,i=1,n+1 
t=1 
Do 30,j=1.n+1 
If(i.ne.j)then 
Š T=(t*(x0-xQ))/(xG)-x0)) 
Else 
Endif 
30 Continue 
Px=px+t*(f(i) 
20 Continue 
Write(*,*)xp 
Stop 
End 


3-4 حذف جاوس: طريقة التعویض الخلفي 


قبل الدخول في التطبیق البرمجي c‏ آری من الأهمية بمکان طرح المفهوم النظري 
hl‏ 4 جاوس أو طريقة ( التعویض الخلفي) . 


الأساسيات النظرية 


تقوم طريقة جاوس على تحویل المعدلات إلى نظام المصفوفات و التعامل معها 
على هذا الاساس كما الاتي : 


А1х1 + а2х2 + a3x3 + و‎ anxn = cl 


B1x1 + b2x2 + b3x3 + ......... bnxn= c2 








XI 


g ka‏ هة فضف AS‏ المجاهیل, و المراد ایجاد قیمها $ بینما ا 








(augmented matrix‏ عن 












all 212 .... aln cl 
b21 2,2 ...... b2n 2 









à; العمود الثاني و هذا‎ a E الصف‎ ¢ p 


ا ЮР‏ = هي ج > عل a‏ 1 
صفراء و ذلك باتباع З‏ 
F2,1-b2,1/a1,1‏ 

b2,1- b2,1 - f2,1*a1,1 





02,2 =b2,2 - 2 
b2,3-b3,2 — f1,2*a1,3 


الذي قمنا به ۱ هو الاتي : 





2 a MIA p j E at. ه‎ l ١ 13 اھ‎ ZX vi 5 ۲ 
A لف )1 د یه‎ XQ a ы 72 مے الب‎ ١ |) کہ‎ 4 «t 
l 1 J! بة خوارزمية | لامج با رم تعريف‎ | 





а(50,50) b (50, 50) f(50) x (50) عرف المصفوفات‎ . 


1 
2 
readn 3 
4 


do . 
from I = 1 ton read elements of matrix a 
m=n*(n-1) 
do .6 
from 1-2 tom 
from k=1 to 1-1 
if (a (kk) 40 yes f(Lik)-a(Lk) / a(k,k) 
do .7 
from j = 1 ton +1 
a(I,j)=a(1,j) - f(Lk)*a(k,j) 


no 


from I =I - 1 toI- 1 
from j=1 to n+1 
bi j= ad (ز,‎ 
a(l,j)=ad+1,j) 
a(l-1,j)-b(l,j) 
do .9 
from I = п to 1 step – 1 
from j-I-1 ton 
sl=a(I,n +1) 
x(i) = (s1 — $2)/ а(Т,Г) 
$2 =0 
do .10 


write x(i) 








(n+1 








المصفوفة ابتداء 








FORTRAN language ( البرنامج بلغة(الفورتران‎ 








Dimension a(50,50) ,b (50, 50 ( , £(50) , x(50) 





Write(*,*)' enter numbers of your elements (п)! 


Read (*,*)n 





Do 10,i=1,n 
Do 20 , j=1 п+1 

Write(*,*)' ас) =" 
Read(*,*)a(Lj) 

20 continue 

10 continue 
do 40,k=1,n-1 
if(a(k,k).ne.0)then 
do 40,j=1.n+1 
a(i.j)-a(Lj)-f(Lk)*a(k,j) 
40 continue 
else 


do 50,L=i-1,i-1 


do 60,j-1,n-1 

БИ, са) 

а(1,))=а(1+1,)) 

a(l+1,j)=b(Lj) 

60 continue 
50 continue 

do 70,i=n,1,-1 

do 80,j=i+1,n 

s2=s2 + a(Lj)*G) 

80 continue 

sl=a(I,n+1) 

x(i)=(s1 — s2)/a(L,i) 

70 continue 

do 90,i=1,n 

90 continue 

stop 


end 





هذا الكود الذي يقوم Lay‏ سلف و الآن لنبدأ في شرح الكو 













يقوم بحل المء ادلات ت لیت تتساوی = عدد المعدلات م on‏ عدد عم B‏ 


ad, j )‏ -(ز,00 
а(1,)=а(1+ 1, j)‏ 


- 





C-a 
A-b 
C-b 





الخلف 2 كالآتي: 







الأخير إلى الصف الأول بسالب خطوة واحدة 
d‏ سق ара‏ الذي رتبته الصف الاكير 
we‏ الأخبر си T ñ‏ اصل المعادلة إنما 





kal TA يي‎ j ; 
a = wk : 





bu 1 ارية من اجل‎ ew» و جملة‎ s2 | LA اله‎ | à i | К 5 7 | T 






FIT 


FITTING METHOD الملائمة والانكفاء‎ l-5 


عندما تعطی بیانات غير دقيقة le gi‏ ما أي تقريبية کالقیاسات المعملية 
والهندسية فان هذه البیانات ترتبط بأخطاء مصدرها القیاس أو الانسان أو . 
لذلك نلائم المنحنی الناتج من هذه البیانات حيث أن هذه البیانات نتبع معادلة معينة 
ولکن Айал‏ المنحنی FISS‏ من شکل ولکن المنحنی المفترش Ады‏ من 
jou a‏ البیانات نفسها هل تملل المسئلة دالة تربيعية أو أسية أو ... وهکذا وبعد 
تحدید الدالة يجب أن نجد آجود Best Fit LL‏ وهي تعني أن نجعل الأخطاء 
آقل مایمکن فنأخذ مجموع مربعات الأخطاء 


S- уг 
1=0 


бг = P.(x)-vi 
رط‎ (X) = ао+аух +....... +ашх" =) а(х) 


وهذه الطريقة هي طريقة ملائمة المنحنیات باستخدام المربعات الصغری 
۳ 05 
وتکون dš S‏ ما يمكن عندما эр?‏ 


وسندخل الآن في مثال عملي على ذلك 

الجدول أسفله يعطي بيانات عن مدى استهلاك الماء بإحدى البلدان وذلك 

ببلايين الجالونات في اليوم 

- استعمل الانكفاء الأسى لاستهلاك الماء بدلالة الزمن 

- استعمل ما حصلت عليه لحساب استهلاك الماء في سنة 1975 وقارن بما 
كان متوقعا وهو 449.7 





























1970 1960 1950 1940 1930 السنة 
411.2 322.9 87 136.43 110.5 | الاستهلاك 
(جدول (l-5‏ 


وسنطبق ما سبق على المعادلة التالية (الاسية) 
Y=ab*‏ 
ولكي تصبح خطية نأخذ لوغارثم الطرفین 
Ln (у) = In ( аб")‏ 
Ln(y) = In(a) + x In(b)‏ 
Z=A+Bx‏ 
B=In(b) =‏ و A=In(a)‏ و Z=In(y)‏ 


s- у`5° _Y(A+Bx-9 
i=0 i0 


05 x 
—=2) A-Bx- = 
A 2 +Bx-Z)()=0 


95 n 
—=2) (A4Bx- =() 
B 3 +Bx-Z)(x) 
المعادلات العمودية‎ 
nA ^ B): X; — У; Zi 


AY X i+ ву хі = 3 ZiXi 





z -In(y) وبارجاع‎ 
nA + ву X; = У In(V;) 
AY x+ BY xí = Y xiln(yi) 











n У xi 
A = sau) xi” — xD unes (1) 
) xi У хі? 
Ууд ES 
= ГА = [У хі" mo- У xi У xiln(yDY 
) xiln(yi) y x 
Аз 2) 
п ino) 
B- ГА = Га У; xiln(y)- У myi У xiy A ....... (3) 
У’ ) xiln(vi) 





а=е^ , b-eP 
b за بالعودة للمعادلة الرئيسية و التعویض بقیم‎ 
































)4( ل و ی ما Y=a‏ 
وطبقا للبيانات المعطاة والتي تبين (сла‏ استهلاك الماء باحدی البلدان ببلایین 
الجالونات. 

1970 1960 1950 1940 1930 السنة 
411.2 322.9 202.7 136.43 110.5 الاستهلاك 
(جدول 2-5( 


ولو آردنا ایجاد الحل یدویا فنکمل الجدول ونختصر السنوات كالتالي: 


























x IV Ln(y | X° X*In(y) X^*In(y) 
30 110.5 4.7 900 141.15 4234.51 
40 136.43 | 4.9 1600 196.63 7865.30 
50 202.7 5.3 2500 265.59 13279.32 
60 322.9 5.7 3600 346.64 20798.43 




















70 411.2 | 6.0 4900 | 421.34 29493.49 
































250 | 1183.73 | 26.6 | 13500 | 1371.35 75671.05 








(جدول 3-5( 


بالتعويض عن ad‏ المجاميع 

A = 5%13500- (250)? = 0 
A= (26.6*13500-250*1371.35)/5000 =3.25 
В = (5*1371.35-250*26.6)/5000 =0.041 

а=е?`?° = 25.79 
b=e"""4 = 1.04 

Y=ab" 
Y=25.79*(1.04)* 
у (75) =25.79* (1.04) 8 


هذه الأرقام العشرية ليست دقيقة بل مقربة لرقم أو رقمين ولکن النتيجة النهائية 
ملف ده من ينات Nae Magus‏ 
الجدول بعد الملائمة 





30 40 Bu 60 70 























104.36 |147.94| 209.72 297.31 421.46 





(جدول 4-5( 
قائمة بالمتغیرات المستخدمة في الخوارزمية والبرنامج 


JS متغیر‎ : 5 

iN‏ عدد النقاط المدخلة 

Xi, yi‏ :4„ ٭٭(السنة) و 7(الاستھلاك) المعطاة 

ج : АА‏ المراد Ам)‏ قيمة الاستهلاك عندها 

۷ : قيمة الاستهلاك النانج عند x‏ 

161۴2: قيمة المحدد العام الناتج من المعادلات السابقة...... المعادلة )1( 


X =Š £ مجمو‎ : 5131112 














(sumx=sumx+xi) хі احسب مجموع‎ .6( 





(sumxy=sumxy+xi*In(yi)) xi*In(yi) مجموع‎ Gus! .64 
1=1+1 احسب‎ .65 

66. هل I<n‏ 13 كان نعم ارجع إلى 4 

67. احسب 448 delta‏ من )1( a‏ من المعادلة ) 

TR .70 











المخطط الانسيابي 
Flow Chart‏ 









sumx=sumx+x 
sumy=sumy+y 


sumxx-sumxx-4x* 


sumxy=sumxy+x*In(y 


T-T1 


<< 
© 
— 





#include<stdio.h> 

#include<math.h> 

void main )( 

{ 

double a,b, sumy, y, sumxy, delta, yi[5]; 








printf ("Enter Points Number: 
scanf ("%d", &n) ; 

printf ("Enter Values of x and у"); 
sumx-0; sumy-0; sumxy-0; sumxx-0; 








for (1=0;і<п;і++) 

{ 

fscanf ("%d", &xi[i]); 
scanf ("%1f", &yi[i]); 





vilil-log(vilil); 

INSURE 

i[il; 

Е х+х1 [1] *xi [1]; 

вуне i #y [1]: 

} 

delta = n*sumxx-sumx*sumx ; 

а=ехр( (sumy*sumxx-sumx*sumxy)/delta); 
b=exp( (n*sumxy-sumx*sumy)/delta ); 









printf("Enter value of x"); 
scanf ("%d", &x) ; 

y=a*pow (b, х); 

P writ ng 


= $1f , b= %1f\n\n",a,b); 
п\п"); 
("x 1"); for (i=0; i<n; i++) printf )" 5058", xi [i] 23 
$ t£ [ИУ 2"); سیت‎ VELEN 14у ртїл®Ё ("45,214 
",а хром (b, xi[i])) 
printf("\n\n y = 35. 21۶ * $5.21f^$d = sifinw',a,b,x,v); 
printf ("When x-$d then y=%5.31f",x,y) ز‎ 


7 

















#include<stdio.h> 

#include<math.h> 

void main )( 

{ 

int i,n,sumx,x, sumxx, х1 [5]; 

double a,b,sumy,y,sumxy, delta, уі [5]; 

FILE *stream; 

// input values 

stream = fopen("fitread.FIL", "r"); 

fscanf (stream, "3d", єп); 

(1=0;1<п;1++) fscanf (stream, "%d",&xi[i]); 

for (1=0;і<п;1і++) fscanf (stream, "512", &у1 [1]); 

fclose (stream); 

// Calculation 

sumx=0; sumy=0; sumxy=0; sumxx=0; 

for (i=0; i<n; i++) {yi[i]= =log (yi [i]); 
sumx=sumx+xi[i]; 
sumy=sumy+yi[i]; 











x-sumxx4xi[i]*xi[i]; 
sumxy-sumxytxi[i]*yi[i]; 
) 
delta - n*sumxx-sumx*sumx ; 
-exp( (sumy*sumxx-sumx*sumxy) /delta); 
хр( (n*sumxy-sumx*sumy)/delta ); 
printf("Enter value of x "); 
scanf("$d",&x); 
y-a*pow (b, x); 
// writing 
stream = fopen("fitwrite.FIL", "wt"); 
fprintf (stream, "a= $1f , b= %1f£\n\n",a,b); 
fprir са " У = а * b ^ x\n\n"); 
2 5 Wy P E $ 
;i<n; of fprintf (stream 
fp: : (stream, "\пу :"); 
for (i=0; i<n; i++) i ی ارو‎ черер "$5.21f ",a*pow(b,xi[i] )); 
fprintf (stream, "\n\n у = %5.21Е ٭‎ $5.21f^$d = 
",a,b,x,y); 
fprintf (stream, "When x-$d then y-$5.31f",x,y); 
fclose (stream); 


) 

























الادخا لات 
Enter Points Number: 5 |‏ 
Enter Value of x 30 40 50 60 70‏ 
Enter Value of y 110.5 136.43 202.7 322.9 411.2‏ 
Enter value of x 75‏ 





النتائج 
a= 36.633592 , b‏ 


1.035513 
yz=za*b*x 


x: 30 40 50 60 70 
y :104.36 147.94 209.72 297.31 421.46 


= 36.63 * 1.04^75 = 501.804287 


у= 





| x=75 then у=501.804 










оа] illas‏ في تلك البلد 
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1 iul 
النقاط المعطاة مع منحنی الملائمة‎ 
(3- 5 (شکل‎ 


استهلاك 


ta 


فی الختام لا بسعنا الا أن نسأل 
الله je‏ وجل تحقیق الفائدة من هذا 
الکتاب ونتمنی أنه لقي استحسانکم . 
ویسرنا أن نتلقی آراءکم حول فحوی 
الكتاب على عنوان البرید الالكتروني : 
tkne@tkne.net‏ 


ولکم جزیل الشکر 





وادارة موقع النقنية 


